Why is the bulk resistivity of topological insulators so small? 



Brian Skinner, Tianran Chen, and B. I. Shklovskii 

Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

(Dated: October 12, 2012) 

As-grown topological insulators (TIs) are typically heavily-doped n-type crystals. Compensation 
by acceptors is used to move the Fermi level to the middle of the band gap, but even then TIs have a 
frustratingly small bulk resistivity. We show that this small resistivity is the result of band bending 
by poorly screened fluctuations in the random Coulomb potential. Using numerical simulations of 
a completely compensated TI, we find that the bulk resistivity has an activation energy of just 0.15 
times the band gap, in good agreement with experimental data. At lower temperatures activated 
transport crosses over to variable range hopping with a relatively large localization length. 



The three-dimensional (3D) topological insulator 
(TI) [IHS] has gapless surface states that are expected to 
exhibit a range of interesting quantum phenomena [SJ [7] . 
While a number of 3D TIs have been identified, most of 
these are poor insulators in the bulk, so that the prop- 
erties of the surface are obscured in transport measure- 
ments. For this reason achieving a bulk-insulating state 
remains an active topic of research [8]? -416], 

Typically, as-grown TI crystals are heavily doped n- 
type semiconductors, and correspondingly they exhibit 
metallic conduction. In order to make them insulating 
these TIs are compensated by acceptors. With increas- 
ing compensation K = Na/No, where No and Na are 
the concentrations of monovalent donors and acceptors, 
repsectively, the Fermi level shifts from the conduction 
band to inside the band gap and then at K > 1 into the 
valence band. When compensation of donors is complete, 
K = 1, the Fermi level is in the middle of the gap and 
the most insulating state of the TI is reached. The hope 
is that at K = 1 the bulk resistivity p should obey the 
activation law, 

p = Po exp(A//c B T), (1) 

with an activation energy A that is equal to half the band 
gap Eg. Here, po is a constant and koT is the thermal 
energy. Since typically E g ~ 0.3 eV, this expectation 
would imply a well-insulating bulk at room temperatures 
and below. 

The typical experimental situation at K — 1, however, 
is frustrating [15] . In the range of temperatures between 
100 and 300 K the resistivity is activated, but with an ac- 
tivation energy that is three times smaller than expected, 
A — 50 meV. At T < 100 K the activated transport is 
replaced by variable range hopping (VRH), character- 
ized by p oc e~xp[(To/T) x ] with x < 1, and the resistivity 
grows even more slowly with decreasing T . In Ref. [15] 
the authors show that Mott VRH (x — 1/4) provides a 
reasonable fit to their data at 50 K < T < 100 K. Def- 
inite characterization of the temperature exponent x is 
difficult, however, due to the relatively narrow window 
of temperature and to variations between samples. At 
T < 50 K the resistivity saturates due to the contribu- 
tion of the surface states. 
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FIG. 1. Energy diagram of a completely compensated TI with 
band gap E g . The upper and the lower straight lines (E c and 
E v ) indicate the unperturbed positions of the bottom of the 
conduction band and the ceiling of the valence band; the mid- 
dle line (p) corresponds to the Fermi level. Meandering lines 
represent the band edges, which are modulated by the fluctu- 
ating potential of charged impurities; R g is the characteristic 
size of these potential fluctuations. The percolation levels for 
electrons, E e , and holes, Eh, are shown by dashed lines; the 
activation energy A corresponds to the difference E e — p (or 
p — Eh)- Puddles occupied by carriers are shaded. Shallow 
impurity levels are not shown because they merge with the 
band edges. 

In this paper we suggest an explanation for the un- 
expectedly small bulk resistivity of TIs. We assume 
that both donors and acceptors are shallow and we use 
the theory of completely compensated semiconductors 
(CCS) [13 IE]. This theory is based on the idea that 
near K = 1, when almost all donors and acceptors are 
charged, random fluctuations in the local concentration 
of impurities result in large fluctuations of charge. The 
resulting Coulomb potential is poorly screened because 
of the vanishing average concentration n = No — Na of 
screening electrons. Huge fluctuations in the random po- 
tential bend the conduction and valence bands edges and 
in some places bring them to the Fermi level, thereby cre- 
ating electron and hole puddles that non-linearly screen 
the random potential. Thus, the amplitude of fluctua- 
tions is limited only by the semiconductor gap E g . As 
a result the ground state of a CCS, shown in Fig. [T] 
is similar to a network of p-n junctions [171 118) . The 
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characteristic size of these p-n junctions, also called the 
nonlinear screening radius, is given by 
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where k is the dielectric constant, e is the electron charge, 
and N = N D = N A . For N — 10 19 cm" 3 and k = 20, 



R g ~ 70 nm 3> N^ 1 ^ 3 « 4.6 nm, so that we deal with a 
very long range potential. As a result, the resistivity can 
be dramatically different from the expectation outlined 
above, which assumed flat bands. First, at relatively high 
temperatures conduction is due to electrons and holes be- 
ing activated from the Fermi level to their corresponding 
classical percolation levels (classical mobility edges), E e 
and Eh, in the conduction and the valence bands. These 
may be substantially closer to the Fermi level \x than 
Eg/2, but so far the resulting value of A has not been 
studied theoretically. Second, at sufficiently low tem- 
peratures electrons and holes can hop (tunnel) between 
distant puddles, so that variable range hopping replaces 
activated transport. In the low temperature limit p(T) 
should obey the Efros-Shklovskii (ES) law of VRH [TS], 



p exp (Tes/T) 



1/2 



(3) 



where T E g = Ce 2 //t£, £ is the localization length of states 
with energy close to the Fermi level, and C is a numerical 
coefficient. So far the magnitude of Tes and the nature 
of the crossover between activated and VRH conduction 
have not been studied. 

In this paper, motivated by the TI resistivity puzzle, 
we return to CCS and model numerically the K = 1 case. 
For moderately large T we find that A = 0.15T g . For 
a TI with E g = 0.3 eV this implies A = 45 meV, in 
agreement with observed values [15] . We also calculate 
the single-particle density of states (DOS) of impurity 
states, and we find that the DOS has a Coulomb gap at 
the Fermi level 19J. We show from our simulation that 
the resistivity is described by Eq. ([3| at low tempera- 
tures and crosses over to Eq. Q at higher T. We present 
a crude estimate of the localization length £ which sug- 
gests that Tes ~ 900 K and that the crossover between 
activation and ES VRH occurs at T ~ 40 K. Together our 
results for the activated and VRH resistivity establish a 
universal upper limit for the resistivity p(T) that one can 
achieve for a 3D TI compensated by shallow inpurities. 

In order to model the CCS numerically, we simulate a 
cube filled by an equal number of randomly positioned 
donors and acceptors (20000 of each). We numerate all 
donors and acceptors by the index i and we define rii = 
0, 1 as the number of electrons residing at impurity i and 
the variable fa to discriminate between donors (/j = I) 
and acceptors (/, = — 1). The resulting Hamiltonian is 
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where qi = (/,* + l)/2 — m is the net charge of site i 
and all energies are defined relative to the middle of the 
band gap. The first term of Eq. Q contains the ener- 
gies of donor and acceptor sites, which for the case of 
shallow impurities is very close to ±E g /2. The second 
term is the total interaction energy of charged impuri- 
ties. For two impurities at a distance r » a^, where 
as is the Bohr radius of impurity states, one can de- 
scribe their interaction using the normal Coulomb law 
V(r) = e 2 /nr. For example, an empty donor shifts the 
energy of an electron at a distant filled donor by an 
amount — e 2 /nr. On the other hand, for a pair of im- 
purities with separation r < a^, quantum mechanical 
averaging over the electron wavefunction becomes impor- 
tant (such close impurity pairs are common in heavily 
doped semiconductors, where as > N^ 1 / 3 ). A pair of 
very close donors, for example, cannot create an electron 
state deeper than that of the helium-like ion with bind- 
ing energy 2e 2 /kclb- In order to capture this quantum 
phenomenon in an approximate way, we use the classical 
Hamiltonian of Eq. Q with a truncated Coulomb poten- 
tial V(r) — e 2 /n(r 2 + a 2 ^) 1 / 2 . The result of this trunca- 
tion is to eliminate the unphysically deep electron states 
that would result from very compact impurity pairs with 
an unmodified I/r interaction. We will show below that 
our results are mostly insensitive to the details of this 
truncation. Note that Eq. Q docs not include the ki- 
netic energy of electrons and holes in the conduction and 
valence bands and, therefore, aims only at describing the 
low temperature (fc^T <C E g ) physics of CCS. 

In all results below we use dimensionless units for r, 
ag, £, H, E gi and kgT, measuring all distances in units 
of N^ 1 / 3 and all energies in units of e 2 N 1 / 3 /n. Thus, Eq. 
Q can be understood as dimensionless, with E g ^> I and 
V(r) = {r 2 +a 2 B )- 1 / 2 . For a TI with E g = 0.3 eV, k = 20 
and N = 10 19 cm -3 , the unit of energy e 2 N 1 / 3 /n w 15 
meV, so that the dimensionless gap E g rj 20. We were 
unable to directly model E g = 20, since in this case the 
very large R g s» 16 leads to large size effects. Instead, 
we present results for the more modest E g = 10, where 
Rg s» 4 and size effects are negligible, and for E g = 
15, where R g « 9 and size effects can be treated using 
extrapolation. Unless otherwise stated, results below use 
as = 2 and are averaged over 100 random initializations 
of the donor and acceptor positions. 

In our simulation, we first search for the set of electron 
occupation numbers {rii} that minimizes H. We start by 
assuming that all donors are empty (n; = 0, Qi = 1) and 
that all acceptors are filled (rt^ = 1, % = —1). These 
charged donors and acceptors create a random Coulomb 
potential whose magnitude exceeds E g . We then sequen- 
tially choose pairs consisting of one filled site and one 
empty site and attempt to transfer an electron from the 
filled site to the empty site. If the proposed move lowers 
the total system energy H , it is accepted, otherwise it is 
rejected. To describe the change in H resulting from such 
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FIG. 2. (Color online) Dimensionless single-electron DOS 
g*(e) = g(e) /[2N/ (e 2 N 1 ^ 3 / 7c)] for a completely-compensated 
semiconductor with as — 2 and £" g = 10. The inset shows the 
DOS near the Fermi level e = (upper curve, blue). For com- 
parison, the quadratic Coulomb gap g{e) = (3/7r)e 2 is shown 
by the dashed line The lower (magenta) line shows 

separately the DOS of rare filled donors and empty acceptors. 

a transfer it is convenient to introduce the single-electron 
energy state, £j, at a given impurity v. 

The process of transferring electrons concludes when all 
pairs i, j with rij = 1 and rij = satisfy the ES stability 
criterion: 

Ej -e i -V(r ij )>0. (6) 

This final arrangement of electrons can be called a 
pseudo-ground state, since higher stability criteria of 
the ground state (involving multiple simultaneous elec- 
tron transfers) are not checked. Such pseudo-ground 
states are known to accurately describe the properties 
of the real ground state at all but extremely small ener- 
gies nana m]. 

Once the pseudo-ground state {rii} is known, the DOS 
g{e) is calculated by making a histogram of the single- 
electron energies {£i}. The result is shown in Fig. [2] 
with the DOS in units of 2N/{e 2 N 1 / 3 /K), so that the 
total area is equal to unity. Occupied and empty states 
are separated by the Fermi level at e — 0. The nearly 
constant DOS between — E g to E g reflects a practically 
uniform distribution of the random potential from —E g /2 
to +E g /2. Near the Fermi level one sees the Coulomb gap 
that is a universal result of the ES stability criterion [19] . 

Once the energies {ei} are calculated, we evaluate the 
resistivity using the approach of the Miller-Abrahams 
resistor network [TSJ. Namely, each pair of impurities 
i,j is said to be connected by the resistance Rij = 
Rq exp[2rtj / £ +Sij / }:bT], where the activation energy 



is denned [TSJ as follows: 

etj = { ( 7 ) 
{ max[|£ 4 | , \sj\] , Ej£i > 0. 

The resistivity of the system as a whole is found using a 
percolation approach. Specifically, we find the minimum 
value R c such that if all resistances Rij with Rij < R c 
are left intact, while others are eliminated (replaced with 
R = oo), then there exists a percolation pathway con- 
necting opposite faces of the simulation volume. The 
system resistivity p{T) is taken to be proportional to R c , 
which captures the exponential term while details of the 
prefactor are ignored [TS]. 

In Fig. [3] we plot the computed resistivity as a func- 
tion of temperature, using the dimensionless logarithm of 
the resistance (lnp)* = (£/2)1ii(R c /Rq) and the dimen- 
sionless temperature T* — 2fcsT/£. These notations are 
introduced to exclude any explicit dependence on £. Fig. 
[j^a) shows {lap)* versus (T*) -1 / 2 over the huge range of 
temperatures 0.03 < T* < 200. One can see that at low 
temperatures T* < 0.3 the resistivity is well described 
by the ES law, Eq. with C w 4.4. The higher tem- 
perature range 1 < T* < 200 is plotted separately as 
a function of 1/T* in Fig. [3^b). Here we find two acti- 
vated regimes of hopping conductivity. At extremely high 
temperatures T* > 50 we see the large activation energy 
E a ~ 0.75E g while in the intermediate range 1 < T* < 10 
we see an activation energy A = {0.15 ± 0.01) E g . We re- 
peated this analysis for the larger band gap E g = 15 
using systems of 10000, 20000 and 30000 donors and by 
extrapolating to infinite size we find A = (0.15±0.02)_E g . 
These results for A remain unchanged, within our statis- 
tical uncertainty, if we use cig — 1 instead of as = 2. 

It should be noted that the large activation energy 
E a ~ 0.75E g observed at T* > 50 does not have any 
physical meaning for a real CCS, since at such large 
temperatures the conduction is not due to hopping but 
rather to free, "hot" carriers far from the conduction and 
valence band edges. Nonetheless, for our model Hamil- 
tonian this result is consistent with established theories 
which say that at such large temperatures E a = (£«), 
where (...) denotes averaging over all pairs i,j (see Ch. 8 
of Ref.HSJ. 

On the other hand, the second activation energy A = 
0.15-B g makes full physics sense and should be seen in 
experiment. At T E g electrons optimize their conduc- 
tivity by hopping among impurities that are energetically 
close to the Fermi level. The activation energy A can be 
understood as the resulting percolation level for hopping 
between nearest-neighboring sites. In other words, if elec- 
trons are activated only to those sites with |e| < e p , then 
precisely at e v > A = 0.15E g there exists an infinite 
conduction pathway for electrons comprised of hops of 
length ~ A^" 1 / 3 or shorter. In a heavily doped semicon- 
ductor this energy is equivalent to the activation energy 
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FIG. 3. (Color online) The temperature dependence of the 
resistivity for E g — 10 (blue dots). The dimensionless resis- 
tivity (hip)* is plotted in (a) against (T*)^ 1 / 2 to illustrate 
that the resistivity follows the ES law at low temperatures, 
and in (b) against (T*)" 1 to show that the resistivity is acti- 
vated at larger T* , with two distinct activation energies. The 
dashed lines (black) are linear best fits. 

of electrons from the Fermi level to the conduction band 
mobility edge E e . (Of course, holes are activated from 
the Fermi level to their percolation level Eh as well.) For 
a typical TI E g = 0.3 eV, so that we get A = 45 meV, 
in good agreement with typical experimental data |15j . 
(We note, however, that recent experiments on Sn-doped 
Bi 2 Te 2 Se have achieved A ~ 125 meV [IB]. Such large 
activation energies may be associated with deep donor 
impurity levels, which go beyond our model.) 

This activation to the percolation level persists un- 
til much smaller temperatures, where A becomes pro- 
hibitively large compared to the thermal energy. At 
such small T* conduction proceeds by VRH among elec- 
tron/hole puddles at the Fermi level and the resistivity 
is given by Eq. 

One can interpret the relatively small numerical fac- 
tor 0.15 above by recalling that in a typical 3D contin- 
uous random potential, ~ 17% of space has a potential 
smaller than the percolation level [18] . As we demon- 
strated above the energy of the conduction band bottom 
is roughly uniformly distributed in the interval (0,E g ). 
This means that the percolation level E e should be close 
to 0.17 Eg and makes our result A = 0.15E g quite rea- 
sonable. 

So far we have emphasized results that do not explicitly 
depend on the localization length £. In fact, knowledge of 
£ is necessary to predict Tes and the transition tempera- 



ture T t between Eq. ([!]) and Eq. ^ in real temperature 
units. (According to Fig. [3^,, the transition happens at 
T* « 1/2, or T t ~ £/4). We argue now that in a TI £ is 
quite large, leading to a prominent role for VRH. To see 
this, consider that if an electron with energy close to the 
Fermi level is assumed to tunnel from one electron puddle 
to another distant puddle along the straight line connect- 
ing them, then the tunneling path passes through regions 
where the conduction band bottom is quite high above 
Fermi level. This implies a small tunneling amplitude, 
or £ <C as- In fact, however, a tunneling electron can 
use the same geometrical path as a classical percolating 
electron with energy A above the Fermi level. In order to 
roughly estimate £, we assume that along such a classical 
percolation path the tunneling barriers V are uniformly 
distributed in the range < V < A and we neglect the 
curvature of this path. Integrating the action along this 
path then gives £ ~ ^/(toA) 1 / 2 = flB\/e 2 /aflA. For a 
TI with E g = 20 and a B = 2 this gives £ ~ 0.8. This 
crude estimate leads to T ES — 900 K and T t ~ 40 K, 
which is similar in magnitude to the experimentally ob- 
served Tt ~ 100 K where the resistivity crosses over from 
activated to VRH behavior IT-T. 

We note that if one plots our result for (bap)* against 
(T*)-V4 in trie relatively narrow crossover range 50 K 
< T < 100 K, one gets a mostly straight line, as seen in 
Ref. [TU However, our results suggest that at low tem- 
peratures the bulk resistivity follows the ES law of VRH 
with temperature exponent x = 1/2, which should be- 
come apparent if the bulk resistivity can be probed to 
very low temperature. Such measurements are presum- 
ably possible in samples that are much thicker than those 
studied in Ref. US (~ 100 //m). For such thick samples 
conduction through the bulk of the TI crystal dominates 
over the surface transport until much smaller tempera- 
tures. 

To conclude, we have studied numerically the bulk re- 
sistivity of a TI crystal with band gap E g as a CCS. We 
find that at high temperatures k B T > 0.03E g the re- 
sistivity is activated with relatively small activation en- 
ergy 0.l5E g , in agreement with experimental data [15j . 
At lower temperatures the resistivity crosses over to 
ES VRH, with an estimated characteristic temperature 
T ES ~ 900 K. Thus, Eq. ^ with T ES - 900 K at T < 40 
K and Eq. Q with A = 0.15E g at T > 40 K give the 
upper limit for resistivity that one can achieve for a heav- 
ily doped and completely compensated TI with shallow 
impurities. 
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